Magnetization plateaux induced by a coupling to the lattice 
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We investigate a simple model of a frustrated spin- 1/2 Heisenberg chain coupled to adiabatic 
phonons under an external magnetic field. Using field theoretic methods complemented by exten- 
sive Density Matrix Renormalisation Group techniques generalized to include self-consistent lattice 
distortions, we show that magnetization plateaux at non-trivial rational values of the magnetization 
can be stabilized by the lattice coupling. We suggest that such a scenario could be relevant for some 
low dimensional frustrated spin-Peierls compounds. 

PACS numbers: 75.10 Jm, 75.10 Pq, 75.60 Ej 



The field of quantum spin chains offers a wonderful 
playground for botli theorists and experimentalists to in- 
vestigate a variety of exotic phases cooperatively induced 
by frustration and magnetic field. The so-called zig-zag 
chain with nearest neighbor (NN) and next nearest neigh- 
bor (NNN) Heisenberg couplings Ji and J2 is a funda- 
mental and simple model of a quantum (S=l/2) spin sys- 
tem exhibiting a quantum (i.e. T — Q) phase transition 
(at zero magnetic field) between a quasi-ordered antifer- 
romagnetic phase (so-called Tomonaga-Luttinger liquid 
(TL)) and a spontaneously dimerized gapped phase 0. 
The related quantum critical point is known accurately to 
be located around (J2/Ji)crit = 0.2411 ^J. Interesting 
enough, SrCu02 0] and copper germanate (CuGeOs) 
are fairly good experimental realizations of the zigzag 
chain in the uniform and dimerized phases respectively. 

The underlying richness of the zig-zag chain physics is 
also manifest under an external magnetic field The 
magnetic phase diagram shows, besides the previously 
discussed dimerized phase (at sufficiently low field) and 
several types of TL phases including a chiral phase 
(with spontaneously broken parity), a new phase which 
exhibits (i) a spontaneous breaking of the lattice sym- 
metry of period q — i and (ii) a magnetization plateau 
at 1/3 of the full moment, M = 1/3 (we normalize 
the magnetization M as being 1 at saturation). Note 
that both features (i) and (ii) are expected simultane- 
ouslyiTom the quantization condition qS{l — M)) inte- 
ger |3- Note that the 1/3 plateau state is only stable 
in the range 0.56 < J2/J1 < 1.25 Recent Density 
Matrix Renormalization Group (DMRG) computations 
suggest also that this state supports fractional magneti- 
zation Sz = ±1/3 domain- wall type excitations 

Although its large variety of different phases the mag- 
netic phase of the simple zig-zag chain model does not 



show other plateau phases besides the 1/3 plateau state. 
In this Letter, supported by both analytical and numer- 
ical calculations, we argue that a moderate lattice cou- 
pling can generate an extremely rich magnetic phase di- 
agram with a zoo of new M = p/q (rational) plateau 
states. Experimentally, the lattice coupling is known to 
be crucial in spin-Peierls materials like GuGeOa '4]. It 
has also been proposed to be responsible for a sponta- 
neous tetramerization in the spin-1/2 LiV205 chain 
compound 11]. A cooperative effect of the magnetic field 
and the coupling to an adiabatic lattice was shown to 
produce in 2-leg spin ladders long-range modulated struc- 
tures I2J for several rational values of the magnetization 
M. Although the quantization condition suggests that 
the modulated states of Ref. 01 could give rise to mag- 
netization plateaux, a theoretical investigation of lattice- 
induced plateau phases in quantum spin systems has not 
been carried out so far [13. In this Letter such an in- 
vestigation is performed in the case of the zig-zag chain 
geometry (which can be smoothly connected to the lad- 
der geometry). 

The Hamiltonian of a frustrated spin chain coupled to 
adiabatic phonons in a magnetic field (i?) is written as, 



i i 



(1) 



H is measured in units where gfiB — ^, is the distortion 
of the bond between site i and K the spring constant 
and the first term corresponds to the elastic energy loss. 
Ji sets the energy scale and we fix Ji = 1 in what follows. 

The spin-lattice coupling Ai is dimensionless so that 
the distortions Si are given in units of the lattice spacing. 
Following 1^ we redefine the coupling strengths as Ai = 



2 



Ai{l/ Ky/"^ , although the modulations 5i depend on Ai 
and K separately. 

Let us start with the simplest limit J2 = 0. In that 
case we know that the system dimerizes for any Ai 0, 
opening a spin gap and hence leading to a plateau at 
Af = in the magnetization curve. In the presence of 
an external magnetic field the lattice distortion adapts 
to rest commensurate at any value of M . 

To recover these well known results within bosoniza- 
tion, let us construct the bosonized version of ^ for 
J2 = 0. In the absence of phonons, we can write the low 
energy Hamiltonian for the spin system as 



(2) 

where 4> is the field dual to the scalar field cj) and it is 
defined in terms of its canonical momentum as dx(f> = H. 
The magnetic field effect enters through the Luttinger 
parameter and the Fermi velocity v, which depend 
on the magnetization M. 

We treat the term oc Ai as a perturbation, which in 
the long wavelength limit takes the form 



(3) 



Ai / dxS{x)p{x) 



Here p{x) is the continuum expression of the energy den- 
sity 



= adx(f> + P ■ cos{2kpx 



(4) 



where kp ~ §(1 — M) and a, /3 are M-dependent con- 
stants and the ellipses indicate higher harmonics [l^ . 
The first term simply shifts the magnetic field h — > 
h — aA\bc, where (5c is the fc = Fourier component of 
the modulation, and the second term gives a contribution 
to the energy which reads 



(5) 



It is then straightforward to conclude that the leading 
instability of the lattice deformation that minimizes the 
energy takes the form 



/3^i / dx6{x) : cos{2kFX + V27r(/)) 



S{x) = SQ{M)cos{2kFx) 



(6) 



where So{M = 0) = 6c- This corresponds to the so- 
called fixed modulation which captures the main qual- 
itative features of the model. This statement can be 
verified by computing the lattice modulations in a self- 
consistent way following which solution can be ap- 
proximated by © plus sub-leading higher harmonics 
contributions. The amplitudes of higher harmonics are 
generically smaller than the leading one and will then not 
be considered in the following bosonization analysis [l6| . 



At zero field this modulation produces a total energy gain 
given by 



Ernodm^)}) = KdoiOy - PAiSoiO) J dx : cos(V2^,? 



while for non-zero field we have 



E, 



.od({<5(x)}) = y<5o(Af)2 
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dx : cos(V27r(5 
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(7) 
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If we assume a smooth variation of So{AI) with M [l8| . 
we can conclude that we need a finite magnetic field to 
start increasing M from zero. In that case we have a 
plateau at M = up to a critical field he, after which 
the magnetization jumps to the value Ms such that the 
product —hcMs/2 is of the order of the contribution to 
the energy due to the modulation Emod, Eq- HI)- Due 
to the presence of the relevant term cx cos(\/27r0), the 
system has a spin gap for all magnetizations. However, 
the situation described above (plateau and jump) occurs 
only around AI — and there are no further plateaux in 
the magnetization curve, in accordance with the numeri- 
cal results. The ground state structure above the Af = 
plateau has been studied extensively (see e.g. [lij and 
references therein) and it comes out that a soliton lat- 
tice with a periodicity 2k p starts to develop. The only 
difference found in the magnetization curve between sim- 
ulations with fixed and adaptive modulation (when the 
lattice deformation is determined from minimizing the to- 
tal energy in a self-consistent iterative form) is a change 
in the order of the transition from A/ = 0, that changes 
from first to second order. 

If we add J2 a different situation can occur, and 
in particular non-trivial plateaux can appear in cer- 
tain regions of the parameter space. Let us analyze 
the case of A/ — 1/3 with a modulation of the form 
6{x) = i5o(l/3) cos(^a::). Combining this modulation 
with the second harmonics of the energy density 7 : 
cos{4:kpx + 2\/2TT(j)) : we obtain an interaction energy 
given by 

- Ai(5o(l/3) J dx(^f3: cos(v^(/)) : 

+7 : cos(2\/2^(?!)) . (9) 

To minimize the energy, the second cosine interaction is 
pinned at the minimum of the first one and hence we 
have again a particular situation for Af = 1/3 since the 
second harmonics becomes commensurate only for this 
value of the magnetization. The presence of a plateau 
at 1 /3 depends on the scaling dimension of the second 
cosine interaction, which depends on J2, and from a first 
order analysis one can estimate that it will be relevant for 
values of J2 close to the couplings in CuGeO^, in which 
J2 ~ 0.24- 0.36Ji. 



This is a new generic mechanism for the appearance of 
a plateau due to the spin-phonon coupling. The novelty 
is that the plateau is not produced by the commensura- 
bility of the main (relevant) harmonics (as for the zero 
magnetization case) but it is due to the commensurability 
of the next-to-leading harmonics, whenever it is relevant. 

Note that a plateau at M = 1/3 is present in the Ji — 
J2 chain without phonons, but in that case, the plateau 
mechanism is the usual one (so called classical, since it is 
well visualized in the Ising limit 0) and it is driven by 
the operator : cos(3-\/27r0) : which needs larger values of 
J2/J1 than in the present case to become relevant. The 
present situation is thus much more favorable, making 
it potentially observable in experiments. Moreover, this 
plateau can be present also in the extreme anisotropic 
XY case. 

To study the transition from and to the plateau at 
Af = 1/3, we are in a similar situation as for the M = 
case in the normal chain discussed above from which we 
conclude that we have jumps in M{h) at both ends of the 
M = 1/3 plateau. It would be interesting to analyze the 
formation of a soliton lattice similar to that appearing 
above M = in the present case. We expect that the 
only modification from fixed to adaptive modulations will 
be again in the order of the transition. 

This analysis can be applied to more general situations, 
e. g. for a single XXZ chain where one can also expect 
a 1/3 plateau in the Ising regime. In this case one would 
need a rather big Ising anisotropy A > 10 for the second 
harmonics to be relevant [l^ . 

A similar situation is found for M — 1/2 for the Ji — J2 
case, where the second cosine in (PJ is now replaced by 
: cos(3-\/27r(^) : and is hence less relevant. In the present 
case a first order estimate hints that the 1/2 plateau 
could occur at moderate values of J2. Notice that this 
third harmonics is responsible for the plateau at 1/3 in 
the Ji — J2 case without phonons 0, |20| . 

We now turn to a numerical analysis of the magne- 
tization process of Hamiltonian We have used the 
DMRG method to obtain the ground state energy E{Sz) 
in each subspace of the Sz operator (the z-component 
of the total spin of the chain) on a finite chain of Ns 
sites with open or periodic boundary conditions (OBC 
or PBC). Furthermore, minimizing E — E{Sz) — H Sz 
we have found the magnetization M — as a function 
of the applied magnetic field H . 

To begin with, we assume a phonon field 6i with a 
fixed periodic modulation 5i = 5qcos{t:(1 + M)i), as in 
the previous analytic treatment. In Fig. ^ we show the 
magnetization as a function of H for three different sys- 
tem sizes and open boundary condition. Parameters are 
— 0.4 (for which no plateau is present in the pure Ji- 

J2 chain 0) and A160 — 0.4. A finite size scaling study 
of the critical fields is shown in the inset of this figure. 
The plateaux widths at M — i and M = i extrapole 
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FIG. 1: M{H) for the 36, 60 and 84 sites for J2 = 0.4 and 
ASo = 0.4 in the case of fixed modulation and OBC. The 

inset shows the finite size scaling of the width of the different 
pl„. 
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FIG. 2: M(H) for 36, 48 and 60 sites for .h = 0.5 and 
A\ = 0.8 in the case of an adaptive modulation and PBC. 
The solid line is an extrapolation to the thermodynamical 
limit (color online). 

to finite (although small) values in the thermodynami- 
cal limit, in agreement with the bosonization analysis. 
Let us proceed in a more general way, assuming periodic 
boundaries and minimizing the total energy with respect 
to all non-equivalent lattice coordinates 8i . We use the 
iterative procedure proposed by Feiguin et al [isj and im- 
plemented within a DMRG approach by Schonfeld ct al 
|2l| . The algorithm has been constructed by using an ini- 
tial (periodic) ansatz for the (5^ parameters and obtaining 
a new set of 8i from the adiabatic equation, 

5, = ii < S,.S,+i > (10) 

with the constraint (5^ = 0. The procedure is iter- 
ated until convergence for the energy and the distortions. 
Obtaining the distortion pattern in all Sz subspaces the 
magnetization curve is then generated. In Fig. we 
show M{B) for J2 = 0.5 and Ax = 0.8 . The plateaux 
at M = -i and M — \ and is clearly seen. A finite size 
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FIG. 3; Region of stability of the 1/3 et 1/2 plateaux ob- 
tained from a finite size scaling of the critical fields as a func- 
tion of ^ in the interval [0.5, 0.7] and for Ai = 0.8. Stability 
of other plateau phases is not excluded. 



scaling analyze give 0.1433 for the width of the plateau 
at M = i and 0.1654 for the one at M = i. 

Finally we performed a careful finite size scaling anal- 
ysis of the regions of stability of the two most robust 
plateaux. For that purpose, it is only necessary to con- 
sider the Sz values around magnetizations 1/3 and 1/2. 
Although we have applied the same iterative procedure 
as discussed previously, here we have restricted ourselves 
(at each step) to distortion patterns which fit within the 
expected supercell a procedure which greatly im- 
proves the convergence towards the optimum configura- 
tion. The "phase diagram" representing the region of 
stability of the M = 1/3 and M = 1/2 plateaus with ^ 
is shown in Fig. (Pj). Note that stability of other rational 
plateau phases suggested by the bosonisation approach or 
by the naive fixed modulation calculation (see Fig. are 
not at all excluded. However, such phases which proba- 
bly have quite narrow widths are difficult to identify on 
small clusters. 

In conclusion, we have described a new mecha- 
nism leading to the formation of rational magnetization 
plateau phases. It involves a subtle interplay between 
magnetic frustration and lattice coupling. Our claims 
are supported by both analytical and numerical calcu- 
lations. We suggest that quasi-one dimensional spin- 
Peierls systems, like CuGeOa j4] and Tetrathiafulvalene- 
AuS4C4(CF3)4 jl^, where both phonons and frustration 
play a role, would be the most natural candidates to ob- 
serve such a phenomenon. 

We thank A. Honecker for helpful discussions. 
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